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We cross-correlate cosmic microwave background (CMB) lensing and galaxy weak lensing maps 
using the Planck 2013 and 2015 data and the 154 deg 2 Canada-France-Hawaii Telescope Lensing 
Survey (CFHTLenS). This measurement probes large-scale structure at intermediate redshifts « 0.9, 
between the high- and low-redshift peaks of the CMB and CFHTLenS lensing kernels, respectively. 
Using the noise properties of these data sets and standard Planck 2015 ACDM cosmological pa¬ 
rameters, we forecast a signal-to-noise ratio ss 4.6 for the cross-correlation. We find that the noise 
level of our actual measurement agrees well with this estimate, but the amplitude of the signal lies 
well below the theoretical prediction. The best-fit amplitudes of our measured cross-correlations are 
A 2013 = 0.48 ± 0.26 and A 2015 = 0.44 ± 0.22 using the 2013 and 2015 Planck CMB lensing maps, 
respectively, where A = 1 corresponds to the fiducial Planck 2015 ACDM prediction. Due to the low 
measured amplitude, the detection significance is moderate (« 2er) and the data are in tension with 
the theoretical prediction (~ 2-2.5 a). The tension is reduced somewhat when compared to predic¬ 
tions using WMAP9 parameters, for which we find A 2013 = 0.56 ±0.30 and A 2015 = 0.52 ±0.26. We 
consider various systematic effects, finding that photometric redshift uncertainties, contamination 
by intrinsic alignments, and effects due to the masking of galaxy clusters in the Planck 2015 CMB 
lensing reconstruction are able to help resolve the tension at a significant level (ss 10% each). An 
overall multiplicative bias in the CFHTLenS shear data could also play a role, which can be tested 
with existing data. We close with forecasts for measurements of the CMB lensing - galaxy lens¬ 
ing cross-correlation using ongoing and future weak lensing surveys, which will definitively test the 
significance of the tension in our results with respect to ACDM. 

PACS numbers: 98.80.-k, 98.62.Sb, 98.70.Vc 


I. INTRODUCTION 

Gravitational lensing by large-scale structure is a 
promising cosmological probe. During their cosmic jour¬ 
ney toward us, photons emitted at cosmological distances 
are deflected by the intervening matter. As a result, 
we see a distorted image of the source light distribution. 
Lensing distortions produce non-gaussianity in maps of 
cosmic microwave background (CMB) temperature and 
polarization anisotropies. Lensed galaxies are magni¬ 
fied in brightness and weakly distorted (sheared) from 
their intrinsic shape. Statistical measurements of CMB 
lensing Ml and galaxy weak lensing 015] have been 
achieved recently, and are now a useful tool for precision 
cosmology [6]. 

The cross-correlation of CMB lensing maps with other 
tracers of large-scale structure can provide additional cos¬ 
mological and astrophysical information. For example, 
cross-correlations with galaxy or quasar density maps 
measure the bias of the objects (e.g, IMU), while cross¬ 
correlations with cosmic infrared background (CIB) or 
thermal Sunyaev-Zel’dovich (tSZ) effect maps provide in¬ 
formation on the complex relationship between the dark 
matter and the baryons in different forms over cosmic 
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time (e.g., dusty star-forming galaxies or hot, ionized 

gas) mm- 

Similarly, cross-correlating CMB lensing and galaxy 
weak lensing maps can provide useful cosmological infor¬ 
mation. While CMB lensing and current galaxy lensing 
surveys are most sensitive to matter fluctuations at dif¬ 
ferent redshifts [z « 1-2 and z < 0.5, respectively), their 
cross-power spectrum is sensitive to large-scale structure 
at intermediate redshifts z ss 0.9. Combining the auto- 
and cross-power spectra can thus provide tomographic in¬ 
formation on the growth of structure. Furthermore, the 
cross-power spectrum is immune to nearly all system¬ 
atic effects that can plague measurements of the lensing 
convergence auto-power spectrum (e.g., the point spread 
function (PSF) correction, for galaxy shapes), since the 
CMB and galaxy lensing surveys are completely indepen¬ 
dent measurements. In fact, the CMB lensing - galaxy 
lensing cross-correlation can be used to measure the mul¬ 
tiplicative bias in galaxy lensing shear maps, thus over¬ 
coming an important systematic in cosmic shear analy¬ 
ses [M El- 

Ref. |T5] (H15) reported the first detection of the cross¬ 
correlation of CMB lensing and galaxy lensing with a 
significance of 4.2<r, using CMB lensing maps from At¬ 
acama Cosmology Telescope (ACT) data and galaxy 
lensing maps from the Canada-France-Hawaii Telescope 
Stripe 82 Survey (CS82). They found best-fit ampli¬ 
tudes A = 0.78 ± 0.18 with respect to a fiducial model 
based on Planck 2013 cosmological parameters, and 
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A = 0.92 ± 0.22 for a model based on Wilkinson Mi¬ 
crowave Anisotropy Probe (WMAP) parameters. They 
also noted that uncertainty in the redshift distribution 
of their source galaxies, determined from cross-matched 
COSMOS redshifts for a small subset of the data, could 
cause 10 — 20% changes in the theoretical prediction. 

In this work, we perform a similar analysis using 
CMB lensing maps from the Planck satellite (2013 
and 2015 data releases^ and galaxy weak lensing data 
from the Canada-France-Hawaii Telescope Lensing Sur¬ 
vey (CFHTLen Sfl CFHTLenS has a similar survey size 
and depth as CS82, and the Planck 2015 CMB lensing 
reconstruction noise is comparable to that in the ACT 
lensing reconstruction used in H15 (but with somewhat 
different ^-dependence due to the different resolutions of 
the two experiments). Therefore, we expect our detec¬ 
tion to be of comparable significance to that found in 
H15. Moreover, since the Planck CMB lensing map cov¬ 
ers nearly the full sky, the outlook for cross-correlations of 
these data with ongoing wide-field galaxy lensing surveys 
(e.g., the Dark Energy Survejj^jand Hyper Suprime-Cam 
Survey is promising. We make predictions for these 
surveys, and also compare our cross-correlation results 
between the two Planck data releases. 

This paper is structured as follows. We first introduce 
the lensing formalism in Sec. |TTJ and describe our data 
analysis in Sec. |III[ We then present our results in Sec.|TV| 
and summarize in Sec. El 


II. FORMALISM 

The lensing convergence is a weighted projection of the 
three-dimensional matter overdensity S = 6p/p along the 
line of sight, 


k{ 6) = l dzW (z)6(x(z)0 1 z), 


(1) 


where x( z ) is the comoving distance to redshift 2 and the 
kernel W{z) indicates the lensing strength at redshift z 
for sources with a redshift distribution dn(z)/dz. For a 
flat universe, 


W(z) 


x 





(1 + z ) x( z ) 

H(z) c 

dn(z s ) x(z s ) - x(z) 
dz s x{Zs) 


( 2 ) 


where is the matter density as a fraction of the 
critical density at 2 = 0, H{z) is the Hubble constant 


1 http: //mu. cosmos. esa.int/web/planck 

2 http://www.cfhtlens.org/ 

3 http://www.darkenergysurvey.org/ 
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FIG. 1. The lensing kernels for the CMB (dashed curve) and 
CFHTLenS galaxies (solid curve), normalized to a maximum 
value of unity. 


at redshift z, with a present-day value Hq , c is the 
speed of light, and z s is the source redshift. Note that 
f 0 dz s dn(z s )/dz s = 1. We hereafter denote the galaxy 
lensing kernel computed with the CFHTLenS source red¬ 
shift distribution as IF' tgal (z). For CMB lensing, there 
is only one source plane at the last scattering surface 
2 * = 1100. Using dn{z s )/dz s = 5d(z s — 2 *), where So is 
the Dirac delta function, the CMB lensing kernel can be 
simplified to 


W“(2) = 


(1 + z) X(z) 

0 H{z) c 


X(z*) - x(z) 
X{z *) 


( 3 ) 


The lensing kernels for the CMB and CFHTLenS 
galaxies are shown in Fig. [l] We discuss the CFHTLenS 
source distribution in detail in the next section. The 
mean redshift weighted by the two lensing kernels is 
2 mean = / o °° W Kcmb W Ksa - l zdz // 0 °° W K ™ b W K ^dz « 0.9. 

In the Limber approximation [19], the CMB lensing- 
galaxy lensing cross-correlation is 


C* cmbvi = f — ^Lw K ^ b (z) W K ^(z)P(k,z), 

Jo c x(z)- 

(4) 

where P (k, 2 ) is the matter power spectrum evaluated 
at wavenumber k = £/x(z) at redshift 2 . For our 
fiducial theoretical calculations, we compute Eq. [4] with 
P (fc, 2 ) from the code nicaearl using the nonlinear power 
spectrum from HALOFIT [20j EQ- For a comparison, 
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we also compute theoretical predictions using the halo 
model (e.g., [22][23]) following the methodology described 
in 35s HH (we simply replace the tSZ signal in their ap¬ 
proach with the CFHTLenS galaxy lensing signal). Since 
the halo model is only expected to be accurate to 5-10% 
precision in this context, we use the nicaea+HALOFIT 
approach when comparing our measurements to theory. 
However, the halo model calculation provides intuition 
about the influence of nonlinear power, as it explicitly 
separates the one-halo and two-halo contributions to the 
cross-power spectrum. Finally, we also compute Eq. [4] 
with the linear matter power spectrum from camtj^ for 
an additional comparison. 

The predicted cross-power spectrum is shown in Fig. [2] 
using Planck 2015 cosmological parameters (column 4 
of Table 3 in Ref. [25]). In particular, f l m = 0.3156 
and ag = 0.831, where ag is the rms amplitude of lin¬ 
ear matter density fluctuations at z = 0 on a scale of 
8 h~ l Mpc. Fig.[2]shows that nonlinear contributions are 
non-negligible for i > 100 and are dominant for t > 500. 
Similar results are seen in the halo model comparison, 
where the one-halo term takes over at t ss 600. Note 
that the total power predicted by the halo model is in 
good agreement with the more accurate HALOFIT cal¬ 
culation. 

To demonstrate the cosmological sensitivity of the 
cross-power spectrum, we vary and as by ±5%, and 
show the results in Fig. [3] On most angular scales, the 
cross-power spectrum shows degeneracy between the two 
parameters, where a larger (smaller) S2 m or a 8 simply 
increases (decreases) the overall amplitude. However, on 
very large angular scales (£ < 30), increasing (decreasing) 
fl m decreases (increases) the power. Thus, in principle 
wide-held galaxy lensing surveys covering large sky frac¬ 
tions can break the degeneracy between the parameters. 
Over the range of angular scales considered in this paper, 
the parameters are completely degenerate. 

Later in the paper, we will also compare our measure¬ 
ments to theoretical calculations using the maximum- 
likelihood WMAP9+eCMB+BAO+iLo parameters [26] 
(see their Table 2). In this case, fi m = 0.282 and 
a 8 = 0.817. 

In order to motivate our data analysis below, we con¬ 
sider a simple forecast for the expected signal-to-noise ra¬ 
tio (SNR) of the Planck 2015 CMB lensing - CFHTLenS 
galaxy lensing cross-correlation. We use the fiducial cal¬ 
culation described above with Planck 2015 cosmological 
parameters to compute the signal. We compute the er¬ 
ror bars using the analytic approximation based on the 
auto-power spectra of the CFHTLenS and Planck lensing 
maps (e.g., Eq. 30 of Ref. na) For Planck, we use the 
sum of the CMB lensing signal and noise power spectra 
provided in the 2015 data release, while for CFHTLenS 
we use the measured auto-power spectrum of the con¬ 
vergence maps (thus including both signal and noise as 


6 http://camb.info 



FIG. 2. Theoretical predictions using Planck 2015 cos¬ 
mological parameters. The thin brown solid curve labeled 
“Smith03+Takahashil2” shows our fiducial theoretical calcu¬ 
lation using nicaea+HALOFIT. The other curves show pre¬ 
dictions using the halo model and the linear matter power 
spectrum, as labeled in the figure. 


well). The maps are described in full detail in the next 
section. Adopting the same sky fraction (140deg 2 ) and 
multipole range (40 < i < 2000) as in our analysis below, 
we obtain a predicted SNR ss 4.6. Since this estimate 
is based on the sky-averaged Planck CMB lensing noise 
power spectrum rather than the actual power spectrum 
in the specific CFHTLenS sky patches, the actual SNR 
is expected to differ slightly. However, as noted in Sec. [TJ 
this forecast is comparable to the H15 SNR ss 4.2 ob¬ 
tained using ACT CMB lensing maps and CS82 galaxy 
lensing maps. 


III. DATA ANALYSIS 

In this work, we use CMB lensing maps from the 
Planck satellite data releases in 2013 and 2015, and 
galaxy lensing maps from the CFHTLenS survey. While 
Planck is a full-sky survey, CFHTLenS covers only 
154 deg 2 . Thus, we cut out regions in the CMB lensing 
maps that match the CFHTLenS fields, and construct 
both CMB and galaxy lensing maps in real space with 
the same resolution of 0.16 arcmin 2 per pixel. Fig. [3] 
shows examples of the CMB and galaxy lensing maps for 
the CFHTLenS W1 field. We then combine the masks 
from both surveys and apply them to all data sets. Fi¬ 
nally, we analyze the cross-power spectrum between the 
two surveys, and present our results and null tests in 
Sec. [TV] 
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FIG. 3. Cosmological sensitivity of the cross-correlation. 
We fix the fiducial cosmology at Planck 2015 cosmological 
parameters, and vary and as by ±5%. 


A. Planck data 

We consider the CMB lensing maps produced by the 
Planck collaboration for both the 2015 and 2013 data re¬ 
leases 0E1- Both maps are based on lensing reconstruc¬ 
tions using quadratic estimators [27]. The 2015 map [6] 
is provided as an estimate of the CMB lensing conver¬ 
gence field, reconstructed using the minimum-variance 
combination of all temperature and polarization estima¬ 
tors applied to CMB component-separated maps from 
the SMICA code [28] . The publicly released map is band- 
limited to the multipole range 8 < £ < 2048. We also 
use the associated mask, which removes regions contam¬ 
inated by emission from the Galaxy and point sources, 
leaving 67.3% of the sky. Note that the mean-field bias 
has already been subtracted from the publicly released 
map, and we thus perform no additional such subtrac¬ 
tion in our analysis. 

For a comparison, we also consider the CMB lensing 
map provided in the 2013 data release [3]. In this case, 
the map is provided as an estimate of the CMB lensing 
potential (f >, reconstructed using only the temperature- 
based quadratic estimator applied to the 143 and 217 
GHz Planck 2013 temperature maps. The reconstruction 
noise levels are roughly twice as large in this map as in the 
2015 map [3j. The map is band-limited to the multipole 
range 10 < £ < 2048. We convert the 4> map into a 
convergence map in harmonic space: 

= (5) 

where IZf is the lensing response function provided in 
the 2013 data release. We then transform the resulting 
convergence map to real space in order to extract the 


data in the CFHTLenS regions. 

We combine the associated 2013 lensing mask with the 
2015 mask, although it appears that the 2015 mask is 
stricter and covers essentially all of the 2013 mask, plus 
additional sky regions. In particular, we note that the 
2015 mask covers tSZ clusters, whereas the 2013 mask 
does not. The 2013 reconstruction masks tSZ clusters in 
the 143 GHz channel, but not in the 217 GHz channel 
(where the tSZ signal is null), and thus in the publicly 
released map constructed from a combination of the two 
channels, lensing signal is included at the location of tSZ 
clusters. Since the 2015 reconstruction is based on the 
SMICA map, which combines all Planck channels, tSZ 
clusters are masked prior to the reconstruction in order 
to avoid biases. We test for effects resulting from the 
cluster masking in Sec. [TV] We also note that biases in 
the Planck CMB lensing reconstruction due to tSZ or 
CIB leakage should be small due to Planck’s resolution 
and noise levels [35] 130] , even with no masking, with the 
possible exception of small scales {£ > 1000) in the lens¬ 
ing map (however, the reconstruction noise is large on 
these scales). 

In order to cross-correlate the Planck CMB lensing 
maps with the CFHTLenS convergence maps, we project 
the relevant regions of the CMB lensing maps (and the 
associated masks) onto flat-sky grids in (RA, Dec). This 
procedure uses a cylindrical equal-area projection imple¬ 
mented in the flipper softwar^] which was developed 
by members of the ACT collaboration. The projection 
is performed at high resolution (HEALPix N s i^ e = 8192) 
in order to minimize any resulting artifacts. We verify 
the accuracy of this procedure by calculating the power 
spectra of simulated maps before and after the projection 
(i.e., in the patch on the sphere and in the flat-sky pro¬ 
jection), finding no measurable differences over the range 
of angular scales considered in this paper. 

B. CFHTLenS data 

The CFHTLenS survey is one of the first large galaxy 
lensing datasets that are publicly available (see also COS¬ 
MOS [31]). It consists of four sky patches located far 
from the Galactic plane, Wl, W2, W3, and W4, with a 
total area of 154 deg 2 and a limiting magnitude *ab < 
24.5. The CFHTLenS data analysis pipeline consists of: 

(1) creation of the galaxy catalogue using SExtractor [32]; 

(2) photometric redshift estimation with a Bayesian pho¬ 
tometric redshift code [33]; and (3) galaxy shape mea¬ 
surements with lensHt 0 |34j A summary of the data 
analysis process is given in Appendix C of Ref. [32;. We 
refer the reader to the CFHTLenS official papers men¬ 
tioned above for more technical details. 

The procedure of our galaxy lensing map construction 
can be found in Ref. [35] • In brief, we apply a cut of 


7 http://www.hep.anl.gov/sdas/flipperDocumentation/ 
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FIG. 4. The CMB (left) and galaxy (right) lensing maps for the CFHTLenS W1 field. The galaxy lensing map is smoothed 
with a a g = 1 arcmin Gaussian kernel. No filter has been applied to the CMB lensing map. Data in the white regions 
are masked out due to bright point sources, such as stars in the CHFTLenS map or radio point sources in the Planck CMB 
temperature maps. 


star_flag = 0 (requiring the object to be a galaxy), 
weight w > 0 (with larger w indicating smaller shear mea¬ 
surement uncertainty), and mask < 1 (see Table B2 in 
Ref. [51] for the meaning of mask values). These cuts 
leave 5.3 million galaxies, 140 deg 2 of sky, and an effec¬ 
tive number density of n ga i = 12.5 arcmin -2 , with 


^-gal — 


i(E turf 

n E, wf 


( 6 ) 


where f 1 is the survey sky area excluding the masked 
regions, and i denotes individual galaxies. 

We then reconstruct the convergence map from shear 
measurements using [36] . 


^gal(^) 



7i(^) + 2 72^), (7) 


where /c ga i, 7 i, and 72 are the convergence and shears in 
Fourier space, and Z is the wavevector with components 
(^ 1 ,^ 2 )- Note that in this reconstruction we correct for 
multiplicative and additive biases on the shear, as given 
in Eqs. (4) and (6) in Ref. [35j- Finally, the maps are in¬ 
verse Fourier-transformed into real space, and smoothed 
with a (Jg = 1 arcmin Gaussian window. Ref. [5] identi¬ 
fied 25% of the 172 CFHTLenS pointings, each ss 1 deg 2 
in size, to have PSF residuals. Including these fields can 
bias the auto-correlation function. However, we include 
these regions in this work, as there is no correlation be¬ 
tween the PSF residuals and the signal or noise in the 
Planck CMB lensing maps, and the additional sky area 
is useful for our analysis. Moreover, no significant change 
was seen in the CFHTLenS convergence power spectrum 
(for l < 7000) or peak counts when including these re¬ 
gions in Ref. ('551 . 
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FIG. 5. CFHTLenS galaxy redshift distributions for the 
sum of PDFs of individual galaxies (red thick line) and for 
the best-fit redshift (black thin line). We use the sum of the 
PDFs to model the CFHTLenS dn/dz in our analysis. 


The redshift distributions for the source galaxies are 
shown in Fig. [ 5 ] for both the cumulative sum of the 
redshift probability distribution functions (PDF) of in¬ 
dividual galaxies, and the histogram of the best-fit red- 
shifts. We adopt the former for our analysis. We note 
that we do not apply a redshift cut to the galaxy sam¬ 
ple. Normally, a redshift cut of z < 1.3 is suggested for 
CFHTLenS galaxies, due to the limited number of spec¬ 
troscopic redshift measurements at high-;; and the lack 
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of a near-infrared band. At z « 1.3, the 4000A Balmer 
break leaves the reddest band (z band), resulting in a 
larger photo-z uncertainty 01331. Practically, galaxies 
at z > 1.3 can still have high-quality shape measure¬ 
ments, and their lensing kernel overlaps more with the 
CMB lensing kernel. Including these galaxies thus en¬ 
hances the expected SNR of the cross-correlation signal. 

To estimate the level of uncertainty due to the photo¬ 
metric redshifts, we first compare theoretical models cal¬ 
culated using two different dn/dz, each computed from 
redshifts randomly drawn from the PDF of individual 
galaxies. The resulting theoretical curves are almost 
identical (< 1% difference). We further investigate the 
potential impact from the inclusion of z > 1.3 galaxies, 
which account for « 15% of our total sample. As pointed 
out in H15, due to the strong overlap of such high-z 
galaxies with the CMB lensing kernel, uncertainties in 
their photometric redshifts can lead to non-negligible un¬ 
certainty in the cross-correlation amplitude. Ref. [37] 
compared the summed PDF of CFHTLenS galaxies to 
a matched COSMOS sample, which is measured with 
30 bands and hence can be considered the “true” PDF, 
and found some discrepancies for galaxies with z > 1.3. 
We use the COSMOS data points from Fig. 2 of (37] for 
the highest redshift bin z = (1.30,7.00], and replace the 
PDFs of our z > 1.3 galaxies with the resulting COS¬ 
MOS PDF. The theoretical model computed using the 
COSMOS-corrected PDF is nearly identical to that com¬ 
puted using the full CFHTLenS PDF, with only a slight 
decrease in the overall amplitude (2%). This change is 
highly subdominant to the statistical error in our mea¬ 
surement. 

However, because COSMOS data can also suffer from 
systematic uncertainty in the high redshift tail (see 
Figs. 8 and 9 of [38]), we consider a final, crude test 
in which all z > 1.3 galaxies are manually moved down 
to z = 1.3. Under this extreme scenario, the ampli¬ 
tude of the theoretical model decreases by 20%. How¬ 
ever, Ref. [38] shows that the errors on the high-redshift 
photo-z are approximately symmetric, and thus it is un¬ 
realistic to expect that all such galaxies should be moved 
to lower redshifts. If some were moved to higher red¬ 
shifts, the theoretical amplitude would increase. In the 
absence of a more precise quantifier of these uncertain¬ 
ties, we conclude that systematic uncertainties in our 
cross-correlation results due to photo-z are on the order 
of w 10%, and at most 20%. 


C. Power spectrum and covariance estimation 

We calculate the cross-power spectrum of the CMB 
lensing and galaxy lensing convergence maps in the 
flat-sky approximation using the pipeline developed for 
Ref. [35]• First, to reduce edge effects, we mask the ten 
pixels nearest the edge of each map. We combine this 
mask with the Planck and CFHTLenS masks described 
above, and then smooth the final mask with a ug — 4 ar- 


cmin Gaussian window. We apply the apodized mask to 
the CMB lensing and galaxy lensing maps, and estimate 
the 2D power spectrum as 

( 8 ) 

where the star denotes complex conjugation. Finally, 
we average over pixels in Fourier space with \t\ £ (£ — 
A£/2, l + At/2), for five linearly spaced bins between 
40 < £ < 2000. We correct for the effect of the mask 
using an appropriate / s k y factor (including the apodiza- 
tion), rather than computing and inverting the full mode¬ 
coupling matrix. Results from Ref. [35] and tests with 
simulations indicate that mask-induced effects only im¬ 
pact the power spectrum at i > 7000, whereas we re¬ 
strict our measurement to £ < 2000 here (due to the 
band-limited Planck lensing maps). 

To estimate the covariance matrix, we cross-correlate 
the CFHTLenS galaxy lensing maps with 100 simulated 
Planck CMB lensing maps (for both the 2013 and 2015 
data, separately). We process these simulated maps 
through the same pipeline as the actual Planck lensing 
maps. We then compute the covariance matrix from the 
100 cross-power spectra. The diagonal components of 
the covariance matrix agree to within 10% with the the¬ 
oretical variance estimated using the auto-power spec¬ 
tra of the Planck and CFHTLenS maps (e.g., Eq. 30 of 
Ref. (15|). The off-diagonal components are relatively 
small, < 5% of the diagonal terms. We use the full co- 
variance matrices for all calculations in our results below. 


IV. RESULTS 
A. Measurement 

Fig. |6] shows the cross-power spectra of the Planck 
CMB lensing and CFHTLenS galaxy lensing maps. We 
use cosmological parameters from either the Planck 
2015 [25] or WMAP9 results [25] to calculate the the¬ 
oretical prediction (shown as solid curves). We find the 
best-fit amplitude A with respect to the theoretical pre¬ 
diction by minimizing 

X 2 = E - AC ?)( C J - AC T) (9) 
ij 

where C d is the cross-power spectrum calculated from 
data, C m is the model calculated using Eq. |4j i and j 
denote the multipole bin (five bins for each CFHTLenS 
field), and C _1 is the inverse of the covariance matrix 
d escribed above . The SNR is calculated as SNR = 

VxLn - Xmodel’ where Xnull = X 2( y A = °) and Xmodel 
is the value for the best-fit amplitude A (i.e., minimum 

x 2 ). 

The best-fit amplitudes are shown in Table [I] Using 
the 2013 Planck lensing map, we find Xnull 2013 = 19.1 
and Xmodel 2013 = 15-6 (for either Planck or WMAP pa¬ 
rameters), corresponding to SNR= 1.9. The probability- 
to-exceed (PTE) of the best-fit model is 0.68. Using the 



7 



0 500 1000 1500 2000 

i 



50 

I 

O 


X 




4 K cmb,2015 X K ga 

2 

l- 



- 

0 ' 





+ + Wl 
+ + W2 

2 + + W3 
+ + W4 

- Hanck 

4 - WMAP 

- - best—fit 




■ 


0 500 1000 1500 2000 

l 



FIG. 6 . Cross-power spectra of Planck CMB lensing and 
CFHTLenS galaxy lensing maps. The top panel shows the re¬ 
sult for the 2013 Planck lensing map, while the bottom panel 
shows the 2015 result. The solid curves are the (unsealed, i.e., 
A = 1) theoretical prediction assuming Planck [5SJ or WMAP 
parameters [26]. The best-fit amplitudes with respect to the 
theory curves are A 2013 = 0.48 ±0.26 and A 2015 = 0.44 ±0.22 
using Planck 2015 parameters (shown in dashed curves), and 
A 2013 = 0.56±0.30 and A 2015 = 0.52±0.26 using WMAP9 pa¬ 
rameters. Data points are for individual helds, and errors are 
from the standard deviation of cross-power spectra between 
100 simulated Planck CMB lensing maps and CFHTLenS 
galaxy lensing maps. The boxes represent the inverse-variance 
weighted sum of the four helds. 


FIG. 7. Null tests for the cross-power spectrum. In the upper 
panel, we cross-correlate the CFHTLenS galaxy lensing maps 
and 100 simulated Planck CMB lensing maps. In the lower 
panel, we cross-correlate the Planck CMB lensing maps with 
500 simulated galaxy lensing maps, obtained by randomly 
rotating the CFHTLenS galaxies. Points are for individual 
helds, and errors are the standard deviation of the simulated 
cross-power spectra divided by y/Nsi^ with N s i m = 100 or 
500 for the upper and lower panels, respectively. The boxes 
represent the inverse-variance weighted sum of the four helds. 
We only show results using the Planck 2015 CMB lensing 
maps. The results from the 2013 maps are also consistent 
with zero. 


2015 Planck lensing map, we find Xnuii.2015 = 17-2, and 
Xmodci 2015 = 13.1 (for either Planck or WMAP parame¬ 
ters), corresponding to SNR= 2.0. The PTE of the best- 
fit model is 0.83. In both cases, the model thus provides 
a good fit to the data. 

To estimate constraints on cosmological pa¬ 
rameters, we assume a power-law dependence 


Ci oc (cr 8 ) x ^ ■ Using the theoretical 

model discussed in Sec. [TTJ we find that x ~ 2 in the 
linear regime (£ < few hundred) and x sa 3 in the 
nonlinear regime (£ > 1000)), with a gradual transi¬ 
tion in between. We find that y ~ 1.3 for £ > 200, 
and rapidly decreases to y ~ —0.5 at low-£. These 
power-law dependences are also apparent in Fig. [3j We 
constrain a combination of parameters crs(fl m /0.27) a , 
which parametrizes the degeneracy between and 





































































































A (Planck parameters) 

A (WMAP parameters) 

2013 

2015 

0.48 ±0.26 

0.44 ± 0.22 

0.56 ±0.30 

0.52 ±0.26 


TABLE I. Best-fit amplitudes for the CMB lensing-galaxy 
lensing cross-power spectrum using Planck CMB lensing 
data (2013 and 2015 releases, labeled by the rows) and the 
CFHTLenS galaxy lensing maps. The column labels denote 
whether the amplitude A is measured with respect to a the¬ 
oretical model computed with Planck 2015 cosmological pa¬ 
rameters [25] or WMAP9 parameters [25] . 

flm- For the cross-correlation considered here, we 
find a = 0.41 for the best-constrained combination. 
Assuming a Gaussian likelihood, we obtain a best-fit 
CT 8 (n m /0.27)°- 41 = 0.63lg'i[g. For reference, we also 
list constraints from Planck primordial CMB measure¬ 
ments |2ol [52] and CFHTLenS cosmic shear data m in 
Table [TT| Our constraint remains the same when using 
a = 0.46 for a direct comparison to Planck 2013 and 
CFHTLenS. Our constraint is consistent with that from 
CFHTLenS, but is in « 2er tension with Planck, as seen 
earlier in the best-fit amplitudes presented in Table [l| 
We show two null test results in Fig. [7] where we cross¬ 
correlate (1) CFHTLenS galaxy lensing maps and 100 
simulated Planck CMB lensing maps, and (2) the CMB 
lensing maps and 500 simulated galaxy lensing noise 
maps, obtained by randomly rotating the CFHTLenS 
galaxies. The error bars are divided by y/N s i m with the 
number of simulations N s j m = 100 or 500 for the two 
cases, respectively. The results are consistent with zero, 
with PTE = 0.53 (2013 maps) and 0.11 (2015 maps) for 
test (1), and PTE = 0.61 for test (2) for both Planck 
releases. 


B. Discussion 

The cross-correlation results present some puzzles. 
The SNR of the measurement (« 2) is substantially be¬ 
low the predicted SNR ss 4.6 computed in Sec. [TT] This 
result is entirely due to the low amplitude of the mea¬ 
sured signal with respect to the theoretical prediction. 
The noise properties are as expected — the error bar on 
the measured amplitude for the Planck 2015 lensing - 
CHFTLenS cross-correlation agrees well with the fore¬ 
cast. We find an error of a a = 0.22 (see Table [i]) , while 
the prediction assuming A = 1 is a a = 0.217. 

Thus, the measured amplitude of the CMB lensing - 
galaxy lensing cross-correlation is in some tension with 
theoretical predictions using standard ACDM. The ten¬ 
sion is most significant for the 2015 Planck lensing map, 
as seen in Table [TJ The measured amplitude (A = 0.44 ± 
0 .22) in this case is in tension with the prediction based 
on Planck 2015 cosmological parameters at the 2.5<r level. 
The tension is somewhat less significant for the 2013 lens¬ 
ing map (2 .Oct), due to its higher noise level and some¬ 
what higher preferred amplitude (A = 0.48 ±0.26). Note 



a 

a s (H m /0.27)“ 

Ref. 

Planck 2013 CMB 

0.46 

0-89Lo.o3 

|39 

Planck 2015 CMB 

0.50 

0-90±°; 0 o2 

|25| 

CFHTLenS Cosmic Shear 

0.46 

r\ <77+0.03 
u - 1 1 -0.04 

[40] 

This work 

0.41 

0 63+ 014 

U.UO —o 19 



TABLE II. Cosmological parameter constraints. 


that the decrease in amplitude from the 2013 to 2015 
map (« 10%) is responsible for the fact that the SNR 
hardly improves when using the latter map, despite the 
lower noise (i.e., ss 15% smaller error bar on A). A 
similar amplitude shift from the 2013 to 2015 Planck 
lensing maps is reported in Ref. E3, who use the same 
datasets as in our analysis, but instead cross-correlate the 
Planck CMB lensing maps with the CFHTLenS galaxy 
number density (rather than lensing convergence). They 
use the cross-correlation to infer the linear bias b of the 
CFHTLenS galaxies, finding b = 1.16+g'[g for the 2013 
release (for the 18.0 < iab < 24.0 CFHTLenS galaxy 
sample), but b = 0.821Q']® for the 2015 release, a de¬ 
crease of ~ 30%. The shift we observe is in the same 
direction as that in Ref. EH, but at very low statistical 
significance. Note that cross-correlation with CIB maps 
at 545 GHz in Ref. )6] does not show evidence of a sig¬ 
nificant shift in amplitude from the 2013 Planck CMB 
lensing map to the 2015 map, which suggests that the 
small shift in our results is simply due to noise. 

There are a number of potential reasons for the ten¬ 
sion between our measured cross-correlation amplitude 
and the ACDM prediction based on Planck 2015 cosmo¬ 
logical parameters. One possibility is that the true val¬ 
ues of tjg and O m are somewhat lower than those found 
in the Planck CMB analysis. We note that the ten¬ 
sion between our results and the theory is somewhat re¬ 
duced when comparing to predictions based on WMAP9 
cosmological parameters (see the second column of Ta¬ 
ble [T|. In this context, we refer the reader to the dis¬ 
cussion in Ref. |25| concerning discrepancies between the 
Planck 2015 CMB-determined cosmological parameters 
and those determined from CFHTLenS shear data (par¬ 
ticularly erg and fi m ). It is possible that modeling issues 
(e.g., the nonlinear power spectrum or dn/dz uncertain¬ 
ties) affecting the weak lensing interpretation could be re¬ 
sponsible, although the lowest multipole bin in our mea¬ 
surement in Fig. [6] (where the theory is mostly in the lin¬ 
ear regime) lies clearly below the Planck 2015 theoretical 
prediction. Finally, we note that H15 also found a best- 
fit amplitude for the ACT CMB lensing - CS82 galaxy 
lensing cross-correlation that was slightly low compared 
to predictions based on Planck cosmological parameters, 
though at smaller significance (1.2er) than seen here. 

There are several systematics that could also be re¬ 
sponsible for the observed low amplitude. Photometric 
redshift uncertainties are an obvious suspect, especially 
at high- 2 . We performed three tests in Sec. |III B| to as¬ 
sess the impact of photo-z uncertainties and found them 
likely to be subdominant, but possibly on the order of 
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10%. With presently available data, we are unable to 
fully capture photo-z systematics in galaxy spectral en¬ 
ergy distribution modeling. To accurately quantify such 
uncertainties, one needs observations extending further 
in the near infrared for high-redshift galaxies. Such a 
test is beyond the scope of this work. Using an extreme 
test in which all z > 1.3 galaxies in our data are moved 
to z = 1.3, we find a rough upper bound of 20% on 
photo-z uncertainties in the cross-correlation amplitude. 
Given that the errors on the high-redshift photo -2 are ap¬ 
proximately symmetric }38j , this extreme test is likely an 
overestimate of the effect. Further investigation in this 
area is clearly needed (as noted in H15). We conclude 
that photo -2 errors alone are unlikely to fully explain the 
observed low amplitude of the cross-correlation, but their 
effects are non-negligible (ss 10%). 

A likely physical effect that contributes to the observed 
low amplitude is the intrinsic alignment (IA) of the fore¬ 
ground galaxy shape and the source shape distortion. If 
a foreground galaxy is located between two overdense 
regions, it can be tidally stretched in a direction per¬ 
pendicular to the major axis of the dark matter distribu¬ 
tion. However, the shearing of source light will be aligned 
with the dark matter major axis, and hence the observed 
power spectrum amplitude will be reduced. Refs. mum 
estimate the suppression due to this effect to be ~ 15% 
for CMB lensing - galaxy lensing cross-correlations. To 
fully account for the « 50% difference seen in our re¬ 
sults compared to standard ACDM solely with IAs, a 
very large IA amplitude would be necessary (sa 3 times 
larger than the conservatively expected level). Thus, this 
effect alone is unlikely to fully explain the discrepancy, 
but could be non-negligibl<B 

Another possible contribution to the observed low am¬ 
plitude relates to the mask used in the construction of 
the Planck 2015 CMB lensing map. In the 2015 CMB 
lensing reconstruction, regions where tSZ clusters are lo¬ 
cated are masked prior to the reconstruction. The map 
thus contains no signal at these locations, but because 
these clusters reside in overdense regions where lensing 
signals are expected, the mask could affect our measure¬ 
ment. In contrast, the 2013 Planck CMB lensing recon¬ 
struction contains signal at the location of tSZ clusters, 
because the analysis includes an independent reconstruc¬ 
tion at 217 GHz, where the tSZ signal is null and no clus¬ 
ter masking is required. (The 2015 analysis is performed 
on a frequency-combined SMICA map, and thus cluster 
masking is needed.) In our fiducial analysis above, we 
combined the 2013 and 2015 masks, and thus tSZ clus¬ 
ters are masked. Since the 2015 lensing map simply does 


Recently, Ref. m presented updated calculations of IA contam¬ 
ination for galaxy lensing - CMB lensing cross-correlations, find¬ 
ing that well-constrained low-redshift contributions were consis¬ 
tent with 10-20% contamination, but that unconstrained high- 
redshift contributions could lead to an overall contamination as 
large as 60%. 


not include signal at the location of tSZ clusters, we can¬ 
not use it to study the effect of this masking. However, 
the 2013 map does include such signal (because of the 217 
GHz reconstruction), and thus we can study the effect of 
the tSZ cluster mask by re-running our analysis on the 
2013 map using only the 2013 lensing mask, which does 
not cover tSZ clusters. 

Performing this analysis, we find A = 0.53 ± 0.25 com¬ 
pared to the model based on Planck 2015 parameters, 
and A = 0.62 ± 0.30 compared to the WMAP9 model. 
These amplitudes are « 10% higher than those found us¬ 
ing the 2015 mask which covers tSZ clusters (0.48 ± 0.26 
and 0.56 ± 0.30, respectively — see Table l|). The sky 
fraction in our analysis changes very little between the 
two masks: / s k y ,20i3//sk y ,20i5 - 1 < 0.02. Thus, the in¬ 
creased amplitude is likely due to the inclusion of ad¬ 
ditional lensing signal at the location of the tSZ clus¬ 
ters. However, with only one realization of the sky and 
a relatively noisy measurement, we cannot rule out the 
possibility that the increased amplitude is simply a fluc¬ 
tuation due to including additional data. Testing this 
effect in a dedicated suite of simulations with correlated 
tSZ and lensing signals is needed for a careful assessment. 
Also, note that because current CMB lensing auto-power 
spectrum measurements are almost entirely in the linear 
regime, this effect is likely much smaller there than in the 
cross-correlation with galaxy lensing studied here, which 
receives important nonlinear contributions over most of 
the relevant multipole range (see Fig. [2]). However, this 
effect could be important for the galaxy number density 
- CMB lensing cross-correlation studied in Ref. m- it 
would not explain the difference that they observe be¬ 
tween the 2013 and 2015 Planck lensing maps, because 
they apply the 2015 (and 2013) Planck lensing masks to 
both maps in their analysis. But this effect would bias 
their derived amplitudes low. We defer a careful assess¬ 
ment to future work, but the results above suggest that 
this tSZ cluster mask systematic could explain part of the 
discrepancy of our measured amplitudes with respect to 
ACDM predictions. 

It is also possible, though very unlikely, that the leak¬ 
age of other secondary anisotropies into the CMB lensing 
map could play a role in our results, since these effects are 
correlated with the lensing field (e.g., HUQ3JHSISi!). As 
noted in Sec. m biases in CMB lensing reconstruction 
due to tSZ or CIB leakage are small for an experiment 
with Planck’s resolution and noise levels ED ED , even 
with no masking of clusters or CIB sources. Since the 
2015 Planck reconstruction uses the frequency-cleaned 
SMICA CMB map and further masks the brightest tSZ 
clusters (as described above), such effects are addition¬ 
ally suppressed. Moreover, most of the CIB emission 
comes from higher redshifts than those probed by the 
CFHTLenS lensing kernel, rendering it even less of a 
worry for our analysis. CMB lensing maps can have resid¬ 
ual kinematic SZ (kSZ) signals, as the kSZ effect has the 
same frequency dependence as the primordial CMB fluc¬ 
tuations. However, this leakage vanishes to first order 
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for the kSZ signal, since the line-of-sight velocity of the 
scattering electrons is equally likely to be positive or neg¬ 
ative. Thus, the lowest-order term that could affect our 
results is the kSZ 2 - weak lensing correlation, which is 
highly subdominant compared to the CMB lensing - weak 
lensing correlation that we measure (note that the kSZ 
signal alone is already a second-order effect). The kSZ 2 
leakage into Planck CMB lensing maps was quantified in 
recent simulations performed in Ref. |45j . who found no 
evidence for an impact of the kSZ on cluster mass es¬ 
timation using CMB lensing. Overall, we find it highly 
unlikely that other secondary anisotropies have induced 
significant biases in our results. 


Finally, it is possible that instrumental systematics 
could account for the low amplitude of our measured 
cross-correlation. Measurements of galaxy ellipticities 
are subject to multiplicative and additive biases arising 
from the PSF and other effects that must be carefully 
calibrated. Painstaking analysis using the GREAT and 
SkyMaker simulations is undertaken in Ref. [331 t° per- 
form this calibration for the CFHTLenS data. While no 
significant additive bias is found (and such a bias would 
be very unlikely to cross-correlate with the CMB lens¬ 
ing maps anyhow), a non-trivial multiplicative bias on 
the measured ellipticities, (1 + to) ~ 0.9-0.95, is mea¬ 
sured using the simulations. The bias is larger for low- 
SNR galaxies (i.e., low Wi in the notation of Sec. Ill B), 
which are often high-z galaxies (a fact of particular rele¬ 
vance for our study). Moreover, the uncertainty on this 
multiplicative bias correction is fairly large, with values 
0-85 < (1 + m) < 1.0 consistent with the calibration 
over a wide range of galaxy SNR and photometric red- 
shift (see Fig. 12 in Ref. 05). The multiplicative bias 
propagates directly to the shear and hence the conver¬ 
gence values, which scale as 1/(1 + m). Thus, if the true 
value of (1 + to) is smaller than found in Ref. [53] • the 
derived convergence values will increase, as will the am¬ 
plitude of the convergence auto-statistics and the CMB 
lensing cross-power spectrum studied here. The authors 
of Ref. 01] note the possibility that the galaxy mod¬ 
els considered in their simulations might not be suffi¬ 
cient to capture the true complexity of actual galaxies, 
which could lead to a systematic error in the calibration 
of (1 + to). It is unlikely to be large enough to fully rec¬ 
oncile the discrepancy seen in our results with respect 
to the predictions, but changes in the derived (1 + m) 
values within the allowed range in Ref. 05 could pro¬ 
duce ~ 5-10% changes in the measured cross-correlation 
amplitude. Clearly, this effect also has important impli¬ 
cations for the previously-discussed tension between the 
Planck 2015 CMB-determined cosmological parameters 
and those determined from CFHTLenS shear data [25] . 
especially since (1 + to) enters quadratically in the shear 
two-point statistics. 


Fortunately, this hypothesis can be tested using exist¬ 
ing data, as noted in earlier analyses mm • One can di¬ 
rectly measure ( 1 +to) by computing the cross-correlation 
of ( 1 ) galaxy lensing maps with maps of galaxy number 


density (preferably using a spectroscopic sample), and 
(2) CMB lensing maps with the same galaxy number den¬ 
sity maps. By taking ratios of the measured cross-power 
spectra, one is left with only a factor of (1 + m ) and 
a geometric factor arising from the lensing kernels mi- 
There is also no cosmic variance if one uses the same 
galaxy sample in both measurements. Note that if our 
measured CMB lensing - galaxy lensing cross-correlation 
had higher SNR, we could attempt such a calibration di¬ 
rectly pH (one could also split the galaxy data by the 
shear weight factor and look for a spurious dependence), 
but given the low SNR, a joint approach with galaxy 
number density maps seems more feasible. We leave 
an assessment of this calibration for future work. The 
primary underlying assumption is that the CMB lensing 
measurements are themselves free of a multiplicative bias 
(i.e., that the quadratic estimators have been properly 
normalized). The use of these methods will substantially 
tighten the cosmological constraints from upcoming weak 
lensing surveys mm- 


V. SUMMARY AND OUTLOOK 

Weak lensing of the CMB and galaxies has recently 
emerged as a powerful tool to constrain cosmological pa¬ 
rameters. In this work, we cross-correlate Planck 2013 
and 2015 CMB lensing maps with CFHTLenS galaxy 
lensing maps, and detect a 2 a signal, despite an ex¬ 
pected significance of 4.6er. Our best-fit amplitudes with 
respect to the theoretical predictions are in ss 2-2. 5a 
tension with standard ACDM, with A 2013 = 0.48 ± 0.26 
and A 2015 = 0.44 ± 0.22 using Planck 2015 parameters. 
The tension is reduced if we assume WMAP9 parame¬ 
ters, with A 2013 = 0.56 ± 0.30 and A 2015 = 0.52 ± 0.26. 
A similar discrepancy (but with smaller significance) is 
also found by H15, where ACT CMB lensing maps and 
CS82 galaxy lensing maps are used. We discuss possi¬ 
ble sources of such power suppression, including intrin¬ 
sic alignments (« 15%) and masking of tSZ clusters in 
the CMB lensing reconstruction (< 10%). In addition, 
photometric redshift uncertainties could affect the cross¬ 
correlation at the S 3 10% level. It is possible that other 
systematics not yet accounted for could also play a role, 
such as the impact of nonlinear evolution or baryons on 
the matter power spectrum, or an overall multiplicative 
bias in the CFHTLenS shear calibration. Taken together, 
the combination of all of these systematic effects can per¬ 
haps explain the tension in our results with respect to 
the Planck 2015 ACDM prediction. However, further de¬ 
tailed analysis is needed to understand these effects at 
the required level of precision. 

Due to the limiting size of CFHTLenS survey, less 
than 1% of the available Planck CMB lensing data are 
used in this work. Therefore, future improvement of 
the cross-correlation lies in larger galaxy weak lensing 
surveys, which fortunately are already ongoing. For 
upcoming galaxy weak lensing surveys (the Dark En- 
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77 -gai [arcmin z ] 

/sky 

Expected SNR 

Ref. 

DES 

10 

0.1 

14.4 

46 

HSC 

20 

0.048 

15.5 

47 

Euclid 

35 

0.2 

34.1 

48] 

LSST 

40 

0.25 

39.5 

|49 


TABLE III. Weak lensing survey specifications and SNR fore¬ 
casts. 

ergy Survey (DES), Euclicj^j the Hyper Suprime-Cam 
Survey (HSC), and the Large Synoptic Survey Tele¬ 
scope^ (LSST)), we estimate the SNR of the cross¬ 
correlation with the Planck 2015 CMB lensing data using 
the methodology discussed in Sec. [IlJ For the CMB lens¬ 
ing auto-power spectrum, we use the signal and noise 
power spectra provided in the Planck 2015 release. For 
the galaxy weak lensing auto-power spectra, we consider 
only shape noise in addition to the cosmological signal. 
The weak lensing survey specifications and SNR forecasts 
are shown in Table m The source galaxy redshift dis¬ 
tributions match those used in J21j, except for DES, for 
which we use the redshift distribution given in [50] . The 
predictions are quite promising, with ss 15<r detections 
expected for DES and HSC, and ~ 35-40cr detections 
expected for Euclid and LSST. The predicted SNR for 
the upcoming surveys will be even higher if one considers 
the CMB lensing maps from Advanced ACT m or other 
high-resolution, ground-based CMB experiments. Thus, 


these ongoing and future surveys will precisely measure 
the CMB lensing - galaxy lensing cross-correlation. In 
addition, the comparison of results from multiple inde¬ 
pendent surveys will allow multiplicative shear systemat- 
ics to be overcome. These measurements will definitively 
determine whether the current tension in our results with 
respect to Planck 2015 ACDM parameters is significant 
or simply a statistical fluctuation. 
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